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Abstract 

An assessment of the magnitude of the rearrangement contribution to the Fermi energy and to the 
binding energy per particle is carried out in symmetric nuclear matter by extending the G-matrix 
framework. The restoration of the thermodynamic consistency or, equivalently, the fulfillment of 
the Hugenholtz-Van Hove theorem, is discussed. 
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I. INTRODUCTION 



The physics of nuclear matter beyond the Hartree-Fock (HF) mean field has been dealt 
with a variety of techniques, such as the Bethe-Brueckner-Goldstone hole-line expansion [|l|, 
the variational approach based on correlated functions and the Green's function Monte 
Carlo method ||, exploiting realistic nucleon-nucleon interactions, like the Bonn or the 
Urbana potentials, as input. 

Basically, these studies focused on two observables: The binding energy per particle and 
the saturation density, experimentally extracted from the semi-empirical mass formula and 
from the high energy nuclear electron scattering experiments. In addition, also the com- 
pression modulus, whose value can be inferred from the excitation energy of the breathing 
mode in nuclei, has been quite extensively explored. 

As it is well-known, all these investigations indicate that the theory fails to account for 
the data. A remedy to this shortcoming is presently sought in the introduction of three-body 
forces ||] and/or in a covariant treatment of the nuclear many-body problem ||. 

However, a successful theory of nuclear matter will be required not only to account for the 
data, but to fulfill as well general theorems, in particular the Hugenholtz-Van Hove (HvH) 
one 0. Indeed, a violation of the latter would signal an inconsistency of the theory both at 
the global and at the single-particle level. 

In the following we shall focus on the latter issue, which has received comparatively less 
attention. According to the HvH theorem, whose validity encompasses all the normal Fermi 
systems at zero temperature, at equilibrium the average energy per particle and the Fermi 
energy should coincide. 

The HvH theorem holds because indeed a real (and not complex) energy can be assigned 
to the particles at the Fermi surface, as proved by Luttinger JF]. Actually, the concept of 
single-particle energy is still approximately tenable in the proximity of the Fermi surface, 
but not away from the latter as experimentally verified with exclusive inelastic electron 
scattering (e,e'p) ||. 

Accordingly, in this paper we investigate 

a) the contributions to the single particle energy (in particular at the Fermi surface) 
beyond the mean field, 
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b) how they affect the total binding energy per particle, 

c) whether they help to fulfill the HvH theorem in nuclear matter. 

These items have been considered in the past by various authors. In nuclear matter 
starting from the old paper of Brueckner and Goldman 0] till the more recent work of 



Baldo et al. 10 1 . In finite nuclei, where the issues are even more delicate, the theme has 



been addressed in several investigations. To mention a few we recall those of Faessler et al. 



see, for example, Ref. |TT|) and of Meldner and Shakin |fL2|| . A comprehensive review of the 



topic has been given by Hogdson [|TJ]. 

From the above studies it appears that, in a given theoretical framework, the problem 
of fulfilling the HvH theorem is far from trivial, both in nuclear matter and in finite nuclei. 
This recognition is the basic motivation for the present work. 

II. FORMALISM AND G-MATRIX RESULTS 

We start from the G-matrix expression for the total energy of nuclear matter at zero 
temperature 

k 2 1 

E = J2 7^~ n ( k u M + nE < fcifc2|G|fcifc 2 > n(k 1} k F )n(k 2 , k F ), (1) 

1 felK2 

k F being the Fermi momentum and the G-matrix obeying the equation 

< k\k 2 \G\k\k 2 >=< kik 2 \V\kik 2 > 

_ <y < kik 2 \ V\k 3 k 4 > (1 - w(fc 3 , k F ))(l - w(fc 4 , k F )) < fc 3 fc 4 |G'|fcifc 2 > 
fc 7t 4 e G (k 3 ,k F ) +e G (h,k F ) - e G (kx,k F ) - e G (k 2 ,k F ) 

where the summation includes the spin-isospin degrees of freedom and the matrix elements 
are meant to be antisymmetrized. We then solve Eq. (|2p self- consistently with respect to 
the single-particle energies appearing in the denominator. These, to be defined below, are 
continuous functions of the energy across the Fermi surface. For the bare interaction V we 
use the Bonn potential flTJ . 



Concerning self-consistency, we recall that in Landau's theory of Fermi liquids the energy 
of a single quasi-particle e(k, k F ) is obtained according to the prescription 

SE 



6n(k, k F ) 



e(k,k F ), (3) 



where n(k,k F ) gives the quasi-particle number in the k state. The distribution function 
n(k, hp) fixes the density of the system and depends upon k F , to be viewed as a parameter. 
When Eq. (|3]) is implemented in the G-matrix framework (but ignoring the dependence upon 
n of the G-matrix elements), it yields 

k 2 

e G (k, k F ) = — + Y,< kk 2 \G\kk 2 > n(k 2 , k F ), (4) 

which we view as the energy of a particle of given momentum k. 

In accord with the theory of Brueckner |I5j , which represent the leading term of the hole- 



line expansion and exactly incorporates the two-particle correlations, we solve then Eq. (0) 
seeking for self-consistency on the basis of Eq. ([|). From a diagrammatic point of view this 
amounts to dress the particles and the two-hole lines of a ladder diagram with a first order 
self-energy containing a G-matrix interaction. 

However, it is known from a long time that, when computed at k — k F , Eq. @ violates 
the HvH theorem, in contrast with the HF theory (see the Appendix). This is indeed seen 
to occur in Fig. [I], where we display, as a function of p = 2k F /3n 2 , the results of our 
calculation of the binding energy per particle ([I]) and of the Fermi energy, as obtained from 
(|j), carried out with the distribution n(k,k F ) = 6(k F — k). It appears from the figure 
that our self-consistent G-matrix yields a good binding energy per particle (-16.1 MeV), 
but at the wrong density (p = 0.25 fm~ 3 , i. e. k F = 1.55 fm _1 ); these results are close to 
those recently obtained in Ref. | i~6fl , also in the Brueckner-Hartree-Fock (BHF) scheme, with 



the v 14 Argonne potential and with a continuous auxiliary potential. In addition, a most 
substantial violation of the HvH theorem (of 17.2 MeV, the Fermi energy being e F = —33.3 



MeV) is seen to occur, in quantitative accord with the findings of Ref. |T0 . 

Note that the above quoted e F corresponds to a BHF potential energy of about 80 MeV: 
Hence, the rearrangement corresponds to a correction of = 20% of the latter. 

Thus, in the spirit of Landau's theory, we carry out more thoroughly the variation @ of 
the energy versus the distribution n(k, k F ) writing 

e(k, k F ) = e G (k, k F ) + \- ^ n(h, k F )n{k 2 , k F ) — < kik 2 \G\kik 2 >, (5) 

the last term on the right hand side being often referred to as the rearrangement contribution. 
This in principle should be computed by making the variation on the right hand side of 
Eq. (0) with respect to n(k 3 , k F ) and n(fc 4 , k F ) and also to the n(k, k F ) entering into the single 




FIG. 1: Symmetric nuclear matter binding energy per particle (solid) and Fermi energy (dashed) 
as a function of the density. The lower and upper dashed lines refer to Fermi energies without and 
with rearrangement contribution, respectively, whereas the solid line has no rearrangement. 



particle energies of the denominator. These variations represent corrections to the Pauli 
operator and to the energies of the initial and intermediate states, respectively, resulting 
from the removal of one particle. 

For a discussion of the associated diagrams and for the convergence of the hole-line 



expansion in general we refer the reader to Refs. |T0|, [17], where the impact of the 
hole-hole ladder diagrams and of the off-shell T-matrix in connection with the HvH theorem 
is explored. 

Here we aim instead to compute the rearrangement contribution to the single particle 
energy at (or in the proximity of) the Fermi surface by using the following procedure in 
carrying out the functional derivative, that is 

e(k w k F , k F ) = e G (k, k F ) + - n ( k ^ k F )n{k 2 , k F )—— < fcifc 2 |C7|fcifc 2 > - 

* klk2 OK F op on{K,K F ) 

(6) 

where p = p(k F ) is the density of the system. Of course, if the (large) volume V enclosing 
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the nuclear matter is kept fixed, then 



dp 1 5 1 V r 5n(k u k F ) 1 

2^ n(ki, k F ) = tj^jtt^ / rffci 7-77 r ; = 77- (7) 



5n(k,k F ) V5n(k,k F )j^ v ' y V (27r) 3 7 5n(k,k F ) V 

In the Appendix we shall show that the above prescription for the functional derivative of 
the energy with respect to the distribution, worked out at constant volume, just yields the 
Fermi energy, whereas, when the number of particles is kept fixed, is proportional to the 
pressure. The latter, of course, should vanish for a system at equilibrium. 
To proceed further we now define 



1 d 

A ( k F) = 7777 n (^i> k F)n(k 2 , k F )——(k 1 , fca|G|fci, fe 2 ), (8) 



which, together with Eq. (|7|), allows us to recast (|) as follows: 

e(kxik F ,k F ) = e G (k,k F ) + p-g-A{k F ) 

= e G (k 1 k F )+ f^liip) A{k F ). (9) 

The above equation, where the G-matrix and the rearrangement contributions to the single 
particle (or, better, quasi-particle) energy are neatly separated, yields the correct Fermi 
energy e(k F , k F ) = e F , being rigorously valid at the Fermi surface. 

We should now compute the impact of the rearrangement on the binding energy per 
particle. In this connection, we are aware that the BHF formula 



-tt — 7777^ n(k,k F ) 
N 2iVY 



k 2 

h e G (k, k F ) 

2m 



(10) 



- linking the total and the single particle energies, — no longer holds when the rearrange- 
ment is included. Furthermore, calculations of the nucleon momentum distribution including 
correlations among nucleons beyond the BHF ones indicate that the corrections to the latter 
tend to be more pronounced near the Fermi surface |19| , where the Landau theory applies. 



Hence, inspired by the Landau theory, we heuristically assume Eq. (||) to be applicable as 
well in the proximity of the Fermi surface, where the concept of quasi-particle is tenable. In 
the conclusions the foundation and the limits of this assumption will be addressed. Here we 
write, as an extension of (|T0|), 

f = f + ^ ) fe ln ")" l ^ t £' ,(MF) 
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FIG. 2: The quantity A(kp), defined in the formula (||) of the text versus the density p = 2k F /2n 2 . 
The crosses show the calculated values, whereas the solid line corresponds to the linear fit A = 
150/9+15 MeV fm. We do not display A(kp) for p < 0.05 fm -3 because here its numerical evaluation 
becomes quite inaccurate. 



Eg 



^ + A(k F ) 



' d_ 
dk F 

' _d_ 



In p 
In p 



kp 



dk k l 



(11) 



-i 



1 



kp 



In the above the summation is meant to be carried out in a restricted region near the Fermi 
surface, namely in the range of momenta kqp < k < kp. In the derivation of (0), the 9- 
function distribution is used to be consistent with the G-matrix calculation. Obviously the 
momentum /cqp should be viewed as a parameter: If it will turn out to be close to kp, then 
the omission of the /c-dependence in the rearrangement contribution to the single-particle 
energy (^) should be expected to have not too serious consequences. 

Thus, the expression (|TT|), beyond the standard G-matrix contribution (namely Eq/N), 
explicitly embodies the rearrangement one as well, but the latter, — and this is important, 
— is reduced by the factor [1 — (/cqp//cf) 3 ]/2 with respect to Eq. (|9]). 
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FIG. 3: Symmetric nuclear matter binding energy per particle (solid) and Fermi energy (dashed) 
as a function of the density. The lower and upper solid lines refer to binding energies without and 
with rearrangement contribution, respectively, whereas the dashed line includes rearrangement. 



III. RESULTS FOR THE REARRANGEMENT ENERGY 



In order to compute the rearrangement contribution, the quantity A(kp) is needed: It 
has been numerically evaluated and it is displayed as a function of p in Fig. fj. It appears 
that A(kp) is always positive, growing linearly with p, and substantial. Actually, it turns 
out to be well approximated by the expression A{kp) = (150p + 15) MeV fm. 

Since the above findings on A{kp) have been obtained with a ^-function distribution, 
for consistency we evaluate [<9(ln p)/dkp]~ 1 with the ^-function as well. An elementary 
calculation yields 

Hence, from Eq. @ it follows that the rearrangement contribution to the Fermi energy 
grows as hp, in accord with Ref. [20], where it was conjectured to grow as kp and found to 
vary as k F in a schematic estimate. 

In Fig. m, we display the Fermi energy including the rearrangement, obtained using (|12|). 
versus the density p = 2kp/3Ti 2 . 

We note: 



8 



i) The occurrence, in correspondence of the minimum of the binding energy, of a positive 
rearrangement contribution of about 25 MeV, larger than the value of 17 MeV obtained in 
Ref. [|1(J in a perturbative scheme. Since in Ref. [Ill only the leading diagrams, beyond the 
BHF ones, were kept, one could view the difference as an estimate of the contribution of the 
higher order terms in the hole-line expansion. However, this estimate should not be taken 



at face value, since in [TO the Paris potential was employed, whereas we use the Bonn one; 

ii)notably, our rearrangement contribution stays remarkably constant to the left of our 
saturation density p ~ 0.25 fm~ 3 (say, in the range 0.20 < p < 0.25 fm~ 3 ), in accord with 



the old finding of Thouless |20 



iii) while in Ref. [HJ the obtained rearrangement contribution was just enough to restore 
the thermodynamic consistency, in our case we "overcure" the HvH violation: Indeed, we 
predict the equality between the Fermi energy and the binding energy to occur at p « 0.17 
fm -3 , i. e. not at equilibrium. However, in [|1(| no change was assumed to occur in the 
ground state energy due to rearrangement. 

Here, as discussed above, we schematically estimate the magnitude of this change by 
resorting to formula (|TT]) and searching for a value of the momentum &qp such to restore 
the validity of the HvH theorem. This turns out to be fulfilled when /cqp = 0.09kp, as it 
appears from Fig. ^, where we display, in correspondence to this value, the binding and 
Fermi energies, both including rearrangement, versus the density. 

We see in the figure that the minimum of the binding energy, while reduced, as expected, 
to a value (-13.15 MeV) higher than the BHF prediction, occurs at a saturation density 
Pmin — 0.19 fm -3 (fci? = 1.41 fm -1 ), in closer contact with the experimental value p exp = 0.17 
fm -3 . 

Concerning the rearrangement effect, at kp = 1.41 fm -1 we obtain a positive contribution 
of m 20 MeV for the Fermi energy and of ~ 1.8 MeV for the binding energy per particle. 
These values might be compared with those obtained long ago by Brueckner et al. fl21~| , 
namely about 12 and 1.5 MeV for the Fermi and mean energies, respectively, and with 



the recent ones of Ref. |18[ , namely about 12 and 1 MeV, respectively. It should however 
be remarked that these values have been obtained with simple, as compared to the Bonn 
potential, interactions. 

Finally, the compression modulus predicted in our framework, via a polynomial interpo- 
lation of the binding energy, turns out to be ~ 150 MeV, significantly larger than the one 
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obtained in a pure BHF scheme (« 120 MeV), as expected owing to the rapid Uf dependence 
of the rearrangement term, but still somewhat lower than the experimental value. 



IV. CONCLUSIONS 

In these concluding remarks we seek for some justification of our empirical procedure in 
computing the rearrangement contribution to the binding energy per particle. 

We start by recalling the link (see (|ToD ) between the single-particle and mean energies, 
which holds only for strict mean filed theories, like HF and BHF. Should these schemes be 



valid, then, as it is well-known fL3| , the separation and the single-particle energies would 
coincide, as it is (nearly) true in atoms (Koopman's theorem, see Ref. [|T~3|| ) . 

In nuclei, of course, HF (which respects HvH) is not applicable and BHF, as shown by 
many calculations, while not unrealistic, fails to fulfill the HvH theorem. Indeed, the amount 
of the failure measures the impossibility of describing the system in terms of independent 
constituents. Actually, beyond the mean field framework the only remaining link between 
the mean and single-particle energies is the one expressed by the HvH theorem itself. 

Our simple approach is based on the premise that the contribution to the system's mean 
energy arising from the proximity of the Fermi surface can still be simply related to the 
energy of individual entities. We identify the latter with Landau's quasi-particles and not 
with the BHF particles, since, as previously mentioned, the BHF framework fails to account 
for the appreciable depletion of the single-particle orbits induced by the strong short-range 
repulsion among the nucleons, a depletion occurring mainly, although not only, near the 
Fermi surface. Indeed, it is established that the repulsion much affects both the momentum 
distribution and the rearrangement energy pOf . 

Accordingly, we use the quasi-particle approximation for the propagator G(k, u) in the 
expression 



£ = I lim / f°° *± e *>n 

N p v ^o+ J (2vr) 3 7-oo 2iri 



k l 1 

— + -£*(&, w) 
2m 2 v ; 



TrG(fc,cj), (13) 



which yields the binding energy per particle. The quasi-particle approximation for G(k, u>) 
is best grasped by starting from the canonical spectral representation, 



... , , , SJk,u') r e F . . SJk,uj') 

G(k,u)= / duo' — ^ — — + / duo — ^ — — . (14) 

Je F UJ — UJ 1 + If] J-oo OO — UJ 1 — IT] 
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It amounts to set for the hole spectral function the expression 

S^(k U ) ~ 1 ^{k)\W{k)\ (15) 

where W(k) is the imaginary part of the quasi-particle self-energy and Z(k) the so-called 
quasi-particle strength. 

In our crude and empirical model, to avoid the introduction of further parameters beyond 
/cqp, we have set Z(k) = 1 and W(k) = 0. Of course, the former should be lower than one, 
but not too much, in order not to spoil the concept of quasi-particle, and the latter should 
be small for the same reason. Thus our quasi-particle propagator differs from the BHF one 
only in the location of the pole, which is moved by the rearrangement contribution. 

Finally, we find it gratifying that the only parameter of our approach, namely fcqp, confines 
the quasi-particle existence to a quite narrow domain close to the Fermi surface. Also 
satisfying is our result that the rearrangement affects the Fermi energy an order of magnitude 
more than the mean energy, a finding on which a general consensus exists. 



APPENDIX A 



We comment here on the formula 

8E[n] 



dE dk F bp 



where 



bn(k,k F ) k - k F dk F dp 8n(k, kp)' 
bp 5 



J2 n (h, k F ) 



fci 



(Al) 



(A2) 



Sn(k, k F ) Sn(k, k F ) 

for the functional derivative of the energy with respect to the distribution function. The 
sum over spin and isospin is understood. 

We first consider the volume V constant. In this case 



Hence 



Sp 1 ^ 5n(ki, k F ) 1 

6n(k, k F ) = 6n(k, k F ) V' 



SE[n] dE dk F 1 dE dk F _ dE 

5n(k, k F ) k~k F ~ dk F d(N/V) V ~ dk^~dN ~ dN 



p = €p. 



(A3) 



(A4) 
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In the case of a free Fermi gas, ( |A4j) yields indeed 

5E[n\ 8E 7T 2 1 7T 2 8 (NZky 

~ dk F 2k F V ~ 2k 2 F dk F \V5 2m~ / 



Sn(k, k F ) 



k=kf 



kp 
2m 



FG 



On the other hand, when N is constant (but V varies) 

6N = E x n N , M k > **) = E Sn (k, k F ) = 0, 
^ bn(k,k F ) V 



since 



5N 



^2 n(ki, k F ) = 1. 



5n(k, k F ) 5n(k, k F ) ^ 
Thus, the vanishing of 6N does not imply the vanishing of SN/Sn(k, k F ). Hence 



5E[n] 



5n(k, kj 



dE dk F 1 dE dk F p 5V 



(AS) 



(A6) 



(A7) 



(A8) 



k - k F dk F dp V dk F dp V 8n(k,k F y 
the second term on the right hand side actually not contributing because V and n(k,kp) 
are varying independently. Accordingly 



SE[n] 



5n{k } kj 



k=ki 



dE 1 dk F 1 
dfoNd(l/V) V 



IdE _ P 
p dV p ' 



(A9) 



which vanishes at equilibrium because so does the pressure P. 
Curiously, from (|A9|) for a perfect Fermi gas it follows 



5E\n] 



(A10) 



5n(k, kp) 

showing that a non confined Fermi gas satisfies the HvH theorem, i. e. it reaches equilibrium, 
at zero density. 

On the other hand, for a translationally invariant Fermi system, with a generic interaction 
V(q), one has from ( |AT^ ) 

5E\n] k F d(E/N) 



5n(k } kj 



0. 



k - k F 3 dkp 

In the HF approximation |22j 

\el G + \ P V$) - J *0{k F -k)J **0(k F - k')V(\k - k>\) 



(All) 



E 

N 



(27 

1 1 r 2k F 



1 _ 3 JL _l_ 1 ( q 



22k F 2\2k f 



and it is an easy matter to carry out the derivative in ( |A11| ), getting 



dk } 



2k f 



(A12) 



(A13) 



12 



Now, in HF the Fermi energy reads 



e F F = e F F G + pV(0) - J ^x~Ak F - k')V(\k F - k' 



(2tt)' 



4° + PV(0) 



2k f 



dqq 2 V(q) 1 



(A14) 



(2tt) 2 7 

Hence, by comparing with ( A13 ), one sees that the HvH theorem is fulfilled in the HF 
approximation no matter what the interaction is, providing the latter is independent of the 
density. 
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